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Abstract. General perturbations of a spherical gas bubble in a compressible 
and inviscid fluid with surface tension were proved in [5], in the linearized 
approximation, to decay exponentially, ~ e~'"*, F > 0, as time advances. 
Formal asymptotic and numerical evidence led to the conjecture that F 
4 ^ exp (^-B^^, where < e <C 1 is the Mach number and A and B 

are positive constants. In this paper, we prove this conjecture and calculate A 
and B to leading order in e. 

1. Introduction and outline 

The detailed non-spherical deformations of gas bubbles in a liquid is a problem 
of great physical interest in fundamental and applied physics; see, for example, 
[31 m [S] and references cited therein. We consider the dynamics of a gas bubble in 
a compressible, inviscid and irrotational fluid with surface tension. This physical 
system has an equilibrium state, consisting of: a spherically symmetric gas bubble 
at constant pressure inside the bubble, and a fluid at constant (lower) pressure and 
vanishing velocity, outside the bubble. 

In [5] the linearized dynamics about such equilibria was studied and it was shown 
that general finite energy perturbations of the spherical equilibrium damp out as 
time advances. In particular, the L°° and local energy norms of the perturbation 
tend to zero exponentially, ^ e~'"*, as t — > cx). Formal asymptotic and numerical 
evidence led to the conjecture that 

1) r e) « 2- <^^P 

Here, < e ^ 1 denotes the Mach number, a dimcnsionless ratio of speeds, and 
We denotes the Weber number, a dimensionless measure of surface tension. A and 
B are positive constants determined in [5] by computer simulation. 



(1) 



The asymptotics ([T]) are in marked contrast to the decay rate of perturbations 
which are spherically symmetric: 

(2) r,adiai(e) = 0{e) 

Thus, non-spherical deformations excite shape modes which lose their energy to the 
fluid and radiate sound waves very slowly. 



The Mach number, e, is the ratio of the bubble-wall radial velocity to the sounds speed in 
the fluid, exterior to the bubble. We shall later set We equal to one and consider the asymptotics 
for small e. 
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Our goal in this article is to present a proof of the above conjecture and a calcu- 
lation of A and B to leading order in e. 

A systematic and detailed discussion of the physical problem, the full nonlinear 
compressible equations and the appropriate linearization is presented in [5]. The 
dynamics of the linearized velocity potential, \E'(a;, t), a; G K'^, and the bubble shape 
perturbation, /3(r2,i), e 5^ is governed by the wave-system: 



(3a) 






e^d^"^ - A* 


= 0, 


r = \x\ > 1 


(3b) dr'f 


= dtf3, 


r = 1 


(3c) dt^ 




-^(2 + A5)/3, r^l 


(3d) {p,Yn 


= . 


\m\ < 1 



Equation pap follows from linearization of the continuity and Euler momentum 
equations, exterior to the equilibrium bubble. Equation (j3bp is the linearization of 
the kinematic boundary condition. Equation (j3cp is the linearization of the Laplace- 
Young boundary condition, stating that the pressure jump across the gas - fluid 
interface is proportional to the mean curvature. Finally, (j3dp is the linearization 
of the statement that the origin of coordinates (the center of the bubble) is in 
a frame of reference moving with the bubble center of mass. The equilibrium 
bubble has been normalized to have unit radius. Here, YJ™ = Yf^ifl), E are 
spherical harmonics, which satisfy — AsFj™ = l{l + l)yj™, ^ > 0, \m\ < I, where 
(l^",y,r')L2(s2) = SwSram'. In particular, Yo°{n) = 

As explained in [5] and now outlined, the L°° and local-energy decay, for e small, 
is controlled by a non-self adjoint eigenvalue problem, which arises by seeking time- 
harmonic solutions of (|3a| -(|3d | : 

^{n,t) = e-'^**A(r,i7), i3in,t) = e-'^'px{n) , 

which satisfy an outgoing radiation condition as r — > oo: 
The Scattering Resonance Spectral Problem 

(4a) 

( A + (eA)=' ) ^-A = 0, r>l 
(4b) dr^x = -iX(3x, r = 1 

(4c) -zAv&A = 37 + Wx, Y,°)Y,° - (2 + As)/3a , r = 1 

(4d) ^'a outgoing r — > oo. 

If A is such that ([4]) has a non-trivial solution, then we call A a (deformation) scat- 
tering resonance energy, and (\I'a,/3a) a corresponding scattering resonance mode. 
Due to the radiation condition (j4d|) . the eigenvalue problem Q is non-selfadjoint. 
The eigenvalues of Q all lie in the open lower half complex plane and for each 
e > are uniformly bounded away from the real axis. In particular, for any e > 



For simplicity, wo have set Ca, the cavitation number, appear in [5] equal to one. 
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there is a scattering resonance eigenvalue A*(e) with Im X^,{e) < and such that 
any scattering resonance eigenvalue of A, satisfies 

Ini A < Im A*(e) < ; 

see Theorem 3.2 of [5]. 

The following theorem, a consequence of Theorem 5.1 in [5], relates A*(e) to the 
rate of decay of perturbations about the spherical bubble equilibrium: 

Theorem 1. There exists eq > 0, such that for all < e < the following holds. 
Consider the initial boundary value problem ([3]). Assume the initial conditions: 

(a) (3it,n) = EiyoEimiKil^rmYrm, where 

\m = o)\\ = E E(i + o'^^ i/3r(o)i<oo 

l>0 \m\<l 

(b) *(t,.T), dt'i/{t,x) =0, t = 0, \x\>l, 

Then, there exists a unique solution ^'(r, fl, t), l3{il, t), defined for r > 1, £7 G S'^ , 
which solves the the initial-boundary value problem. 

Furthermore, define r(e) to be the minimum distance of a scattering resonance 
( in the lower half plane ) to the real axis. That is, 

(5) r(e) = I Im A.(e) | 
Then, the solution '^{x,t), /3{^l,t) satisfies the 

Decay Estimate: 

(6) \Pin,t)\ < c = 0)11 e-r(^)*, nes^ 



(7) \^{x,t)\ < 



0, 



Our main result is the following precise asymptotic statement about r(e): 
Theorem 2. For < e < cq, there exist positive constants A,, and B, such that 

(8) r(e) = — exp 



£2 J 

where B w 0.26924 and where A^ has asymptotic expansion 
(9) A, = Ao + O , where Ao « e-^ ^^^ 



Acknowledgement: O.C. and S.T. were supported in part by NSF grant DMS- 
1108794. M.I.W. was supported in part by NSF grant DMS-10-08855. 
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2. Proof of Theorem [2] 

We begin by deriving a characterization of the scattering resonance energies of 
the eigenvalue problem ^ as the zeros of an analytic function in C. 

We first note that outgoing solutions of the three-dimensional Helmholtz equa- 
tion are linear combinations of solutions of the form ft.p'' (r)YJ'"(rj), |m| < I, where 
denotes the outgoing spherical Hankel function of order /. Thus, we seek solu- 
tions of (j4]) of the form: 

where a and b are constants to be determined. This choice of vfi;^ solves the 
Helmholtz equation and satisfies the outgoing radiation condition. To impose the 
boundary conditions at ?■ = 1 we substitute the expressions for and /3\ into 
(j4|) and obtain the following two linear homogeneous equations for the unknown 
constants a and b: 

eA/ip^'(eA) iX \ f \ ^ f ^ 



A^ h^\eX) + (i±M^.AM^)'(.A) 



= 0. 



[ -^Xh\'\eX) -il + 2){l-l)/We j U J " I ' ' ' ^ 2. 
Setting the determinant equal to zero yields 

( 

Finally, multiplying through by and defining 

(11) Q(/,e) = (/ + 2)(;-l)-^ 
yields the following transcendental equation 

(12) zh[^\z) + Q{l,e)h[^'>'{z) ==0, z^eX^O. 
Here h^^"^ denotes spherical Hankel function of the first kind of order /. 



Remark 1. We will analyze the roots z of (|T2| for We ~ 1 since the formula for 
the roots z for We ^ 1 are obtained from those with We = 1 by replacing by 
■jI?^. Furthermore, since X—^, proving Theorem [2] is equivalent to showing that 
the roots of for We = 1 satisfy the property 

A \ B 

(13) r(e) = inf {— Im z\ ~ exp — - 

for < e < eo sufficiently small. 

In the following, we consider the roots in different regimes in / relative to e. 

Definition 3. / = 0{g) means that for any e sufficiently small, there exists a 
constant C independent of e so that \f \ < C\g\. We write f = Os{g), i.e. f is 

strictly of order g, if {j = 0{g) and g = 0{f)j . We write f ^ g if f ~ o{g), that 

*s f / 9 0, as e — > and f ^ g i/ — > oo as e — > 0. 

Outline of the proof: In section [3] equation (|T2|) is rewritten as a coupled set of 
equations for the real and imaginary parts of z = x + iy. In section [3] we show, via 
Lemma |3] and Lemma [5l that the scattering resonance of minimal imaginary part 
does not occur for I = 0{1) or for I ^ e^^. In section[5]we show that the scattering 



LIFETIME OF SHAPE OSCILLATIONS 



5 



resonance of minimal imaginary part does not occur for 1 <C ^ ^ e^^. Therefore, 
the scattering resonance of interest must satisfy I = 0,5(6^^). The detailed analysis 
of this regime and the rigorous approximation of ([S]) is in section |51 Section [7| is an 
appendix containing many of the asymptotic forms of special functions. 



3. Derivation of a system of equation for x, y, where z ^ x + iy 

It turns out that for large / some of the roots z ~ x+iy of (|12p have exponentially 
small y (with also y/x is exponentially small) and therefore the asymptotic analysis 
is more delicate. We next re-express ([T2|) as an equivalent system of equations 
obtained by setting its real and imaginary parts equal to zero. We note that 
(14) 

h^^H^) = + iyi{z), where j/(z) = ^j^Ji+i/^iz) > Vii^) = '\/^-^'+i/2(^)' 

where Ji+i/2{z) and i^;+i/2(2) arc Bcssel functions of order I + 1/2. Substituting 
into ((T^ leads to 

(15) {x + iy)ji{x + iy) + i{x + iy)yi{x + iy) + Qj'i{x + iy) + iQy[{x + iy) = 
Defining the real valued functions 



(16) 



M{x,y) = \[yi{x + iy) + yi{x - iy)] , (x, y) = [yi {x + iy) - yi{x - iy)] 
2 2iy 

^3{x,y) = - [y'lix + iy) + y[{x ~ iy)] , A4{x,y) = —[y'i{x + iy)-y'i{x-iy)], 
2 Ziy 



(17) 

Bi{x,y) = -[ji{x + iy)+ji{x~iy)] , B2{x,y) = -^[jiix + iy) - ji{x - iy)] 
Bz[x,y) = ]^[i[{x + iy) + i'i{x ~ iy)] , Bi{x,y) = -^[j[{x + iy) ~ j[{x ~ iy)], 

and separating out the imaginary and real parts of (1151) . we get 

(18) xAi = -QA^ - y {xBo. + Bi + QB^ - yA^) 

(19) y (^1 + xA2 + QAi) = {xBi + QB3 - y'Sa) 
We rewrite (HH), ([T9|) in the form 

QA3 y{xB2 + Bi + QBi - j/As) 



(20) 



.2^^ _ xBi + QB3 - y^B2 

^ ' ^ Ai+xA2 + QAi 
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4. Location of the roots z of ([12 



As will be seen, the roots with smallest |Im z\ occur for I = 0(e~^), z ~ 0{l); for 
them, it turns out that Imz is exponentially small in e. First, we show that roots 
z in other regimes have larger |Im z\. 

Lemma 4. Ifl = 0{1) as e 0, then any root of hl'^) in the sector arg z £ (— tt, tt) 
satisfies z^ ^ {I + 1)Q, and y = Imz = 0(e^'+^) . 



i (see 



See Theorem 6.2 of [S] for detailed asymptotics. 

Proof. Assume first that wc had z ^ 1. Since in this regime ^ jj'-* 

in the Appendix), it follows that z ^ —iQ = O(e^), a contradiction. Similarly, 
z = Os(l) implies z = O(e^), a contradiction again. Thus z <^ 1. Using the well- 
known expansions of of h[^\z), hf^^ {z) for small z (jST]) and ^ in the Appendix) 
we get 



(22) 



.(1)' 



-Q 



hf\z) 



,y[(f) 
M)Q 



MM 



■ijijz) 



2z^ 



(23) 

implying by iteration 



(2;-i)(/ + i) (; + i)(2Z-3)(2;-3)^ 



+ 0(zV^') 



l + O 



(24) z^ = il + l)Q 



1 



Q 



21-1 (2/-3)(2Z-l)2 



o(qH-\q'+'/^i-'-^/^ 



which is real and positive up to and including 0{Q^) Therefore, x can be computed 
within errors of 0{Q^) = 0(e^') from the series ([24| . Using ()7.ip and the small x ex- 
pansion of Bessel functions (see Appendix), the result follows from straightforward 
calculations. | 



Remark 2. The expression ((24|) holds also in the regime 1 <C / <C e ^ 

Lemma 5. If I ^ ^ 1, then any root of f7^) in the sector argz e (— 7r,7r) 
satisfies z ^ —iQ, implying —Imz = 0(J?e^) ^ e"^. 

Proof. We first claim that there are no roots z <^ I. Indeed, otherwise the expansion 
of h\^\z) and h\^^ (z) in this regime (see (|57)) . (|58)) in the Appendix) would lead 
to a contradiction: z = 0{\/TQ) = 0{l^/^e) ^ L On the other hand, the existence 
of roots z = Os{l) implies, using the asymptotics of h^^^ and (ft.|^'')' (see §7.3- §7.5 
in the Appendix for results in different ranges of , ^ ) in (jl2[) . a contradiction: 

^^^^^^^ = (^y^ = Osile^) ^ 1- Thus, |z| L Using the asymptotic behavior 
of h^i^\z) and its derivative in this regime leads (see section W72\ to z ~ —iQ. | 
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5. Analysis of the case 1 < ; <c e ^ 



Lemma 6. // 1 <C / <C e ^, then any root of II 12]) is given by \24% implying that 
z^OS^/^e)^l^e-~^. 



Proof. Using the asymptotics of h''^\z) and in p^ . the assumption 

z ^ I leads to a contradiction: 2 ^ —iQ = 0{Pe^) ^ /. Assuming z — Os{l) and 
using the asymptotics of and {h^^^y (Appendix, §7.3- §7.5) in ((T^ . also ends 
up in a contradiction: ^^^^^^^ = Os ^yj = Os(^e^) ^ 1- In the regime z = o{l), 

exploiting the asymptotics of /ip''(z), we obtain ((22)) and therefore (p4|) . | 

By Lemma [HI any root of in this case satisfies z = o{l). From (see 
Remark [2]), it follows that the roots are close to the real axis. In this section, we 
estimate y ~ Im z in this regime. 

Lemma 7. Any root z — x + iy of hi 2]) in the regime 1 <^ satisfies 

2'+l r ^ -|/+l/2 



,(1) 



(25) y 



vW+1) /egi/2(/ + 1)1/2 



2e 



2? 



1 - 



Q 



(21-1) 



0{Q'/V 



Proof. Since z ^ Z, we use ([M)) and the expansions of Bessel functions in this 
regime (see Appendix) to conclude 



(26) 



1 pe\2i+i 
Ye \ 2l 



l + O 
1 /eQi/2(; + 1)1/2 



2e 



In the same way. 
Using (|7.ip . one obtains (|25 



= 2 + 



2? 
1 
I 



2;+i 



Q 



l + x 



A. 



(2/-1) 
QA3 



0(QV^' 



/+1/2 



^1 Ai 



= 2 + 0{Q/l) 



Remark 3. Since Q^'^/l^/^ scales like l^'^e < 1, we have \y\ - CI {el^/^f^^^ , 
which is a decreasing function of I in the regime 1 <C ^ <C e~^. Therefore, |y| is 
minimal at the right end, and this minimal value turns out to be larger than the 
values obtained in the next section. 

6. Analysis of the case / = C>s(e~^) 

From (dJ), it follows that when I = Os(e"^), Q = Osie-"^). li z = 0,(Z), but 
|z - l| > /i/^, from the asymptotics and ([73]), it follows that (HH) results in 



(28) 



Hence \y\ 



-iQ\ 



- 1 



|Imz| 3> 1. In the crossover regime, 1 



= 0(/-2/3)^ the 



asymptotics of Bessel functions involve the Airy functions Ai and Bi ([76]) and ([77])) 
whose ratio is Os(l), and thus y cannot be exponentially small. 
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So, we restrict to the sub-region 



< 1 and 1 - 



> r^/^. In this 



regime, since "^jf^ is exponentiahy small in /, by (j7.ip . any root z must have y 
exponentially small. We introduce the scaled variables 

L 



(29) VW+T) = Q = he-^q , a; = -^C , -e' log [^e^y] = V 

From the definition of Q, it follows that (for We =1), 

(30) 



2e^ 

q = l,- — 



L 



From ([201) and (f7T|) . we can replace (|7J|) by 



where 



(31) y 



D 



D = Ai+xA2+QAi 



1 [xAi + gAg + y {xB2 +Bi+ QB4, - 1/^2)] 



Note: the square-bracketed quantity on the previous line vanishes by (|7.ip . With 
the notation in ([291), dSO]) and ^ imply 



(32) 



(33) 



qArt e 



-77/e 



Ai eMi 
= 77 + log 



.{xBi+QB3-y^B2) 



D 



Lemma 8. For any fixed (5 e (0, 1), if G 



C e (^, 1 - ^) and 



V > Vo = 2^* — t'^dt, then there exists eo such that for |e| < eo, F and 

G are smooth functions of (, 77, and e, with 

(34) F{(:,r^;0)=C-h^^^^ , G{Cv;0)=V^2h ^^^^ 



c 



C 



e=0 



7^ Furthermore, in these intervals, the system of equations 



and det^^ 

{F(^,7;;e) = 0, G'(C,77;e) = 0} has a unique smooth solution (C(/*, e), 77(Z*, e)), for 
sufficiently small e, with asymptotic behavior in e determined implicitly from 

(35) Z* = ^==-— 757:77^ + 0(6), 



(36) ?7 = 2/ 



2(1 - c2)3/2 

(t-'-l)i/'dt-e2log[?,C] 



be 



/^(1-2C 



2n 



2(1 -C^) 



2(4 



Proof. Using a standard integral representation of ji ( + 1)C) and ?/; ( 1/1(1 + 1) C) : 
([M)) and ([55)1 . and noting that ^^(Z + 1) = Ue~'^, smoothness of F and G in e fol- 
lows provided Ai and D in (|32|) and ()7.2p are nonzero, which is the case, as shall 
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be seen, for e = 0. Routine calculations, using the asymptotics (|68| -(|73 l) show that 
(37) 

1-2(2 



Ai Vi[x) C 

Therefore, 

(38) 



2/,(l- (2)3/2 



F[C,V]^) = C-l* 7 -777^^ — 7¥r + 0{<^ ), 



C 



2C(1-C2) 



from which the formula for F{C, ry; 0) in follows. Now, setting = in 
immediately leads to psp . For the second part, from the asymptotics of j;, yi in 
this regime, (|68)) - ((73)) . we get 



(39) xBi + QB3 - y^B2 = xji (x) 



CMx) 



xji{x) 



C 2(2(1 - (2) 



(40) 

Therefore, it follows that 



^1±^ _ 1 [xA^ + Q^3 + 2/ {xB2 + Bi + QB4 - 2/A2)] 



yi 



1 2h I T, 



/*(1-2C2) 



2(1 -(2) C 



2(2 



0(.2) 



(41) G((,77;e) =77 + 6^ log ^ -e 



2xBi + QB3-y^B3 



D 



= ri + e log 



hCiiix) 



Vi{x) 



log 



1 + ^-V^ 

(2 



-e^log 



1 , 2/. ^——^ /.(I -2(2) 

2(l-(2)^(2^ ^ 2(2^13^ 



+0(6^) 



From equation (|73p in the Appendix, we obtain 



(42) 



2^=exp 
Vi{x) 



(t-2 _ 1)1/2^^ 



(43) G((,77;e) =77-2;*y^ {r 
2, 1 2;, 



e^logU^.( 



1 + /, 



2(4 



0(6^ 



from which the expression of G((, rj; 0) in follows. The Jacobian of F((, 77; 0), 
G(C,ri;0) with respect to ((,77) is clearly nonzero. The statement of existence 
and uniqueness of (((/*, e), ?/(/*, e)) now follows from the implicit function theorem. 
Furthermore, using (|43l) and dsn in F 0, G immediately implies §6^. | 
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Note 4. We now comment on why the mterval restrictions in Lemma|8]for /*, C and 
r] do not matter for finding the smallest \y\ as a function of l^. First, if < 



then for sufficiently small S, the conclusions of Lemma [7] hold and the corresponding 
\y\ is not as small as implied by Lemma |5] for 1^. in the given interval. On the other 



hand if L > 



for sufficiently small i5, then the conclusions of Lemma[5]hold 



where y = Im z is no longer small. Thus, the restriction e 



is appropriate. Furthermore, for in this interval, there is no need to consider the 
possibility ^ {6,1 — 6) since we have already shown that I = Og , implies 
X — Os (e^^) and argued at the outset that we only need to consider 1 — C ^ l""^^^ . 
When (1 — C^) is small, the asymptotic behavior of is given by a ratio of Airy 

functions, as is the case for 1 — C = 0(^~^/'^) (see ([75)1 and ([77)) ) and the resulting 
|y| is not as small as obtained in ([5^ below, as a consequence of Lemma [51 Also, 
we need not worry about possible solutions for which 77 < 770, since we are seeking 
to maximize 77 (minimizing \y\). 

6.1. MELximizing 77 as a function of and determination of r(e) of ([S|). 
Since Im z ^: y ^ — e^^exp [~^], minimizing \y\ as a function of corresponds 
to maximizing rj. Now, psp implies that Z* increases with Q in any compact subset 
of (i5, 1 — 5) for all sufficiently small e. For now, we consider as a continuous 
variable for the purposes of finding the maximum value of 77. We will show later 
that the maximal value of 77 to the order calculated is the same if Z* takes discrete 
values: ^Hl + l)e^ for I e N. 

We seek a critical value Z* for which di^ 7/ = and df i r] < implying a maximum 
of rj. We note that 



(44) = a,.77 = 2 / {t 



-i)'/'dt-^CiA^-C')'/'--r--i^ 



1 



2(1 -C^) 



2(4 



-e^di, log 



We also note that 



(45) 



0. 



5|) implies 

(1 - c')-V2 



(1 + 2C^)V1 



C(2-C2) 2C(2-C^)2 

Substituting ([35|) and ([JS]) into ([33|), we obtain 
(46) 

1 r .. d-C')'/' 



r2 



2{l~e) 



O(e^) 







1 dt- 



v/W^(2C4 + 5e - 4) 



{2-0 



2C{2 - C 



2Y2 



It is clear that the solution ( of g{C,; e) = has the behavior 



(47) 
where 

(48) 



C = Cm,o + e\m,2 + 0(e") 



Idt 



U,0 



2(4 



g(Cm,o;0) =: go(Cm,o)- 



7n,0) 



-0{e') 



+0(e4)=:.g(C;e) 
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Since 
(49) 

it follows that 
(50) Cm,2 = 



i~(2-C^)- 

C(2 - 



\J~^ Cm, o(2Cm,0 + 5Cm,0 ~ 4) 

2.9;(C™,o)C,o(2-C,o)' 



(2^,0 



5Cm.O 



4) 



4C„,o (2-(2-C^o)2) 



We now seek to find C,m,o which corresponds to the maximal rj. We note from 
that c/() < for C e (o, \/2 - ^2^ =: J and 5^, > for C G (\/2 - v^, . Since we 
are seeking a maximum for 77, we must have 5q < at the corresponding C (noting 
that > 0). Thus the roots of 50(C) = only need to be sought for ( £ J. At the 
right end of the interval J, explicit integration -in terms of elementary functions- 
gives < (~ —0.0678), while go — > +00 when ^ — )- 0. Since go is monotonic in 
J, there exists unique (m,o satisfying go(Cm,o) = 0. The explicit calculation of go 
gives go{0.5S) > and 170(0. 59) < 0, implying that the unique maximum C,„,o is in 
(0.58,0.59)-in fact Cm.o = 0.58134... Using ^ we get Cm,2 « -1.1743... Equation 
P5p also implies that the critical (~ l*,m) that maximizes rj is given by 



(51) U, 



''m,0 



Cm,lCm,o(2 - Cm n) (1 ~ 2C,2„ q) 



0(6^ 



x/l-Cm,02 ■ ^ L (l-C,o)^/' 2(l-C_o)3/2_ 

=: ^m.o + e^Zm,2 + 0(e4) « 0.41535 - 2.40716^ + 0{t^) 
Therefore, it follows that the maximized r\ — r]rn satisfies 
(52) 



= 2U,o f (t-2_l)l/2rft+g2 1 2^^^ 2 / {t'^ - l)^^^dt - 2;,„,oCrn.2 (Cm'o " 



1/2 



log[Zm,oCm,o] + log 



^m,oi^ 2Cm.o) 



2(1 - Co) 



2Cm,0 



=: f?m,o+e^7?m,2+0(e'^) 



0.26924+ 2. 1465e2 + 0(e4) 



Therefore the minimum value oi \y\ is 

(53) |ym| = e^^exp [-e"^?7m,o - Vrn.2] [l + 0(e^] 

« -e-^cxp [-0.26924e-2 - 2.1465] [l + O^e^] 

The proof of Theorem [2] now follows from Remark [TJ 

Remark 5. The calculation ([5^ assumed to be a continuous variable rather 
than discrete, = y/l(T+l)e^ with / G N. This however makes no difference to 
the leading order result in ([53)) . Maximization of rj over / € N results in an optimal 
h that is different from the computed L^,^m by O(e^). Since r/ is a smooth function 
of for — h^„i = O(e^), ??(/*) — ?7(/*,m) = O (/* — Z*,™)^ = 0(e'') and therefore 
discreteness does not affect the leading order asymptotic result ((53l) . 
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7. Appendix: Asymptotics of ji and yi in different regimes 



The results quoted below arc standard and either given in standard references 
such as [T], [5] or follow directly from them. 

The spherical Bessel functions ji{z) = ^/j^Ji+i/2{z) and yi{z) — ■\/5^/+i/2(^) 
have the following integral representations which immediately follows from equa- 
tions (10.9.6), (10.9.7) in [1] (see also http;//dlmf.nist.gov/10.9) using = 1 + 1/2 



(54) Mz) 



(55) yiiz) 



1 

27rz 



1 

27rz 







'i-iy 


+1 Teos 




Jo 



(I + —)t — z sinr 



dT 



cxp 



cxp 



-zsinht — {I + -)t 



dT 



^zsinht + {I + -)t 



dt] 



dt 



which results in the asymptotic representations below by standard Laplace method 
for asymptotics of integrals [2]. The Bessel functions ji{z), yi{z) and the Hankel 
function h[^\z) = ji{z) + iyi{z) satisfy 



(56) 



1 



l{l + l) 



u = 



7.1. The regime z <C Z: We only present the asymptotic regimes relevant to this 
analysis. For z <C / we have ( see equations (10.1.2) and (10.1.3), page 437 of [1] ) 



(57) 
(58) 



(2^ + 1)! 



l + O 



{2l-l)\ 



2'-i(/- l)!z'+i 



O 



2{2l-l) 8(2/- l)(2Z-3) V(' + l)^ 
(Three orders are indeed needed because of cancellations.) For Z ^ z ^ 1, using 
T{x + 1) - 72^6-^x^+1/2 for lar ge X we note the simplification: 



(59) 
(60) 



3i{z) 



2-3/2 

I 



ze \ ' 
21 



1 + 



e /ze 
'h/2 V27 



-i-i 



2{2l-l) 2{2l-l) 8{2l - l)(2l - 3) 
Since the asymptotics is differentiable [3, this implies 



(61) 
(62) 
(63) 



yiiz) 



jM 



-(^ + 1) + 



yiiz) 



{l + l){l + 2)z- 



21-1 (2/-3)(2Z-l)2 
, 2/ + 1 2z2 



O 



21-1 [21 - 2,){2l - If 



z" 
O 



LIFETIME OF SHAPE OSCILLATIONS 

^2 ■ 
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(64) 



ji{z) _ 1 /ze\2'+i 



, z 



yi{z) 2e\2l 
7.2. The regime I, argz e (— 7r,7r). In this case 



(65) 
(66) 
(67) 



h^^\z) ~ - exp 
z 



ji{z) - COS z - 



a + 

I [ Z TT 



(l + l) 



yi{z) - - sin z - 



(^ + 1) 



7.3. The regime Z > 1, ^ 



^ - OsH), 1 - lei > \y\ « 

arg(l - e ("7^/3, 7r/3). From dS]) and jSS]), by tlie saddle point metliod ([2]), 
we obtain 

1 



(68) ji{z) 
where 

(69) u{0 = 

(70) yi{z)^- 



,-1/4 



2VIpTT)e 

-?l/2(? + 1)1/2 {t-^_lf'^dt 



X cxp 



i + !f)+o(r2; 



_ i-y/T^ 

24(1 - ^2)3/2 ^ _ ^2-)l/2 ^ 16 '"^ I 1 + yT"r^ 



1 1 , 



-2 -,\-l/4 



V«(^TT)e 



X cxp 



/l/2(/ + 1)1/2 / [t-2_l]l/2^i 



i_!^ + o(r2) 



(71) 



(72) 



y!(2) 



+ 0(r2) 



hi this regime, we have 
(73) = cxp 

yi{z) 

7.4. The regime Z > 1, f = — 

V'c+i) 

(— 7r/3, 7r/3). hi this case we have 



-2^Ji{i + i) f (t-2 _ lyn-dt 



Os{l), lei-i » |y| « 1, arg(e-i) e 



(74) Ji{z) 



(75) y,(z) 



1 



(i-r')"'/^cos 



y^GTT) / ^(r^F^di - 



(1-e 



sm 



(? + l)7r 
2 

il + l)TT 
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7.5. The regime ^ > 1, ^ = + = Os{l), 1 - ? = 0(^2/3). Here the 

asymptotic behavior is given by 

(76) ,,(,)^^g^!Ai(^V32l/3[l_^f 

(77) ,,(,)^-^^Bi(/2/32V3[l_e]) 
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